Draft version February 1, 2008 

Preprint typeset using lAT^^X style cmulateapj 



WHERE ARE THE HIGH VELOCITY CLOUDS? 

ViCENT QUILIS AND BeN MOORE 
Physics Department, University of Durham, Durham DHl 3LE, UK 
Draft version February 1, 2008 

ABSTRACT 

Recent observations of high velocity clouds (HVCs) have revealed compression fronts and tail shaped 
features of HI suggesting that they are interacting with an external medium. We perform 3-D hydro- 
dynamical simulations of HVCs moving through a diffuse hot gaseous component, investigating the 
behaviour of both extra-galactic dark matter dominated HVCs and nearby pure gas clouds that may be 
accreting onto the Galactic disk via a Galactic fountain. Both scenarios can give rise to similar features 
as observed if the external medium has a density > 10^'' cm^"^. Observations suggest that this may be 
too high for a hot ionised halo or intergalactic gas and supports the Galactic fountain origin, a model 
that can also account for the high fraction of HVCs with tails and asymmetrical morphologies. 

Subject headings: dark matter — galaxies: halos — galaxies: formation — galaxies: evolution 



I. INTRODUCTION 

The nature of the high-velocity clouds (HVCs) seen in 
21 cm HI emission remains controversial since their dis- 
covery over three decades ago (Muller, Oort & Raimond 
1963). Most fundamentally, the distances, hence mass and 
environment of these objects remain poorly constrained, 
largely due to the difficulty of finding background sources 
(against which the clouds would appear in absorption) at 
known, comparable distances. It is possible that the HVCs 
may not have a single common origin, but rather they 
are a mix of relatively nearby (|z| < lOkpc), "Galactic 
fountain" clouds, of material stripped from the Magellanic 
Clouds and other satellites of the Milky Way, and perhaps 
of extra-galactic objects (Blitz et al. 1999; Braun & Bur- 
ton 1999, 2000). This latter possibility is most intriguing 
since the HVCs would provide a link to the numerous dark 
matter substructure halos expected in the hierarchical cold 
dark matter model (Moore et al. 1999, Klypin et al. 1999). 

Previous studies have used the morphological, two phase 
core-envelope structure of the HVCs to infer approximate 
distances (Wolfire et al. 1995). Bland-Hawthorn & Mal- 
oney (1999) have proposed the use of sensitive emission 
measurements, together with a model for the ionizing pho- 
ton flux as a function of position within the Galactic halo 
to infer their distances. A handful of HVCs have upper 
distance limits that are <10 kpc due to the presence of 
background stars (c.f. Gibson et al. 2000). 

Briins et al. (2000, 2001) found that « 20% of HVCs 
exhibit position- velocity gradients, termed head-tail HVCs 
by analogy with the structures of comets. A likely cause 
for their elongated morphology is that they are interact- 
ing with a more diffuse, ionized ambient wind. Further 
evidence for such cloud-halo interactions are evident from 
the compressed leading fronts of many HVCs. In this let- 
ter we use high resolution hydro-dynamical simulations to 
examine the morphologies of HVCs moving through a dif- 
fuse medium. This allows us to infer the efficiency of the 
stripping process and will place limits on the possible envi- 
ronment of the HVCs. A detailed study of the interaction 
between clouds and a diffuse medium can be found in Mur- 
ray et al (1993). 



2. THE STRUCTURE OF THE HVCS 

The typical angular size of HVCs are 9^1° and in 
order to be included in existing catalogues (eg, Briins 
et al. 2000) their HI column density, Nhi, should be 
>10^^ cni^^. For a cloud at distance D, we can esti- 
mate its radius, i?, mean spatial HI density, ric, and mass, 
Mc, as follows. R = 9D ^ 2kpc(6'/l°)(i:>/100kpc) and 
ric = Nhi/2R, i.e. 

"="^""''^'""'(io^^)(^) (to^) 

The cloud mass Mc — fHlNniTrR^mHi , can be written; 

JHi o ^iQW J \1° J VlOOkpcy 

where /hi is the neutral hydrogen fraction. A cloud is 
gravitationally bound when it satisfies the virial relation 
R<T^ ~ GMc- With our model assumptions and CTc ^ 10 
km/s, this translates to 

1012 ^ 2x10^°^ ( ) (^-] 

s2 Jhi s2 V 1019 cm~2 J \l° J \100 kpc ; 

The general conclusion is that the nearby HVCs at dis- 
tances D = 1 — 10 kpc, have masses Mc = 10 — lOOM© 
and must be in pressure equilibrium with the warm ionised 
medium located within and above the galactic disk. HVCs 
at several hundred kpc may also be pressure confined by 
some external medium, or they may be gravitationally 
bound by a dark matter halo. For both type of clouds 
we assume that they consist mainly of warm gas with 
temperature Tc ~ 10^ K, velocity dispersion CTc ~ 10 
km/s, and density Uc- The dark matter dominated HVCs 
are constructed by adopting a dark matter density profile 
p{r) = A/(r2 -I- r^) where A = v'^^^/AttG, Vdr = 15 km/s 
and Tc — I kpc. The HI lies primarily within the constant 
density core of the dark matter halo and we adjust its 
temperature profile to maintain hydro-static equilibrium. 
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A typical cloud may be at a distance of 300 kpc, have a 
radius of w 1 kpc and a total gas mass of 1.6 x IO^Mq. 
This leads to a central column density of 10^" cm~^ which 
is higher than the average quoted by Wakker and van Wo- 
erden (1991) but similar to those observed by Briins et al. 
(2001). 

The diffuse halo material (number density riw) is taken 
to be at a temperature Th ~ 10^ K, comparable to the 
virial temperature for a circular velocity of 200 km/s. The 

wind velocity should be of this same order; this is con- 
firmed by the observed distribution of HVC galactocen- 
tric velocities. Clouds ejected from the disk in a Galac- 
tic fountain would fall back to the disk with a velocity 

« VGMdW^ ~ 200km/s, where Mdisk « 5 x IO^'Mq 
and vertical height ^ w 10 kpc. 

3. THE DENSITY OF DIFFUSE GAS IN GALACTIC HALOS 

The properties of the diffuse outer halo of the Milky 

Way is poorly known. Such a halo is expected to exist 
as left over material that has been unable to cool since 
the formation of the Milky Way. Observations so far pro- 
vide only upper limits on its mean density. Constraints 
come from the measured soft X-ray background (Snowden 
et al. 1997, Benson ct al. 2000), from the dispersion mea- 
sures of pulsars in globular clusters and in the Magellanic 
Clouds, from the kinematics (Moore & Davis 1994) or life- 
time (Murali 2000) of the Magellanic Stream. All of these 
constraints place the density of the ionised Galactic halo 
gas n„(50kpc)<5 x 10-\;m--\ Blitz & Robisliaw (2000) 
derive a halo gas density that is a factor of 2 lower than 
this at a distance of 200 kpc by applying ram pressure 
stripping to the halo dwarf spheroidals. Deep ROSAT ob- 
servations of nearby disk galaxies provide complimentary 
upper limits on a very extended component, but they give 
positive detections of ionised gas close to the disks that has 
a density of n^(lOkpc) » 10~^cm~^ (Bregman & Houck 
1997, Ehle et al. 1998). We caution that direct observa- 
tional limits of the diffuse gas density at 200-500 kpc do 
not exist. 

The fact that the HVCs are moving at 150-300 km/s 
through this medium provides an upper limit to its den- 
sity by considering their terminal velocity (e.f. Benjamin 
& Danly 1997). The equation of motion of a pure gas 
cloud of mass rric, area A and velocity Vc through a diffuse 
medium of density Pw is mdvc/dt = pu:{r)vf.A — mcg{r). 
Assuming that the gas is at terminal velocity will give 
an upper Hmit to Pw{r)- In this case dvc/dt = and 
vt = \/ Nhi ■ (GMjr^ )/ pyj{r), therefore a pure gas cloud 
with column density Nhi = 10^" cm^^, at distances 
of 10 kpc and 100 kpc moving at vt — 300 km/s con- 
strains the local IGM nt„(10kpc) < 3 x 10~^ cm~^ and 
n,„(100kpc) < 3 X 10~^ cm~^ (for pure gas clouds). 

The amount of gas stripped from a dark matter domi- 
nated HVC can be estimated using the analytic approach 
described in Gunn & Gott (1972). The gaseous compo- 
nent will be stripped when the ram pressure of the envi- 
ronment is larger than the gravitational restoration force, 
Pwvfi, = {dcf>jjf^, /dr)a where a is the surface density. This 
is plotted in Figure 1 for our fiducial dark matter dom- 
inated HVC and predicts that it will survive intact for 
halo densities riw < 10~^ cm^"* but would be completely 
stripped if nyj is an order of magnitude larger than this. 



Figure 1 also shows the results for gas clouds embedded 
within a cuspy density profile with p{r) oc r^^. The 
steeper potential can retain gas more easily and may ex- 
plain why tiny dwarf galaxies such as Fornax can hold on 
to gas and form new stars at its center. 
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Fig. 1. — The restoration force at different radii within a dark 
matter dominated HVC. The thin solid curve is for dark matter 
halos with constant density core with radius 1 kpc and the thick 
curve shows the behaviour of a density profile with a central cusp 
p(r) oc r^^. The far right location where the horizontal lines in- 
tersect the restoration force indicates the stripping radius beyond 
which all the gas will be removed. 

4. RESULTS 

The simulations have been tackled as pure hydro- 
dynamical processes. The HI is described by a perfect fluid 
with equation of state p = (7 — l)pe where p is the density, 
e the specific internal energy, and the adiabatic exponent 
7 = 5/3. The pure gas HVCs are constructed to be in pres- 
sure equilibrium with the external hot medium, whereas 
the dark matter dominated clouds have an additional po- 
tential field that maintains their dynamical equilibrium. In 
order to integrate the hydro-dynamic equations describing 
the evolution of HVCs moving through an ionised medium 
we have used the high-resolution shock-capturing method 
described in Quilis, Ibanez and Saez (1996). This code was 
specially designed for an accurate treatment of fluid dy- 
namical processes, and is extremely good at dealing with 
shocks, strong discontinuities, turbulent regions and low 
density regimes. A specially adapted version for the HVC 
simulations has been developed. We employ a fixed grid 
with ratio of sides 2:1:1 and a wind blowing down the long 
axis. The timestep is scit by the; minimum Courant condi- 
tion in the simulation volume which typically sets several 
hundred steps per crossing time. 

We ran a grid of simulations, varying the density and 
temperature of the external medium from 10^'' to 10^® 
cm~^ and 10^ — lO^'K respectively, and cloud velocities be- 
tween 170 km/s and 400 km/s in order to be conservative. 
Most of our simulations were run without allowing radia- 
tive cooling, although a test run with cooling switched on 
showed no difference in the response of the HVC to the 
wind. Our model HVCs are simplistic since they have 
no internal structure, however we do not expect that this 
would greatly influence our conclusions. 
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Fig. 2. — The column density for the dark matter dominated HVC after 300 Myrs of evolution in a wind of velocity 200 km/s, (a) density 
10~^ cm""^ and (b) density 10~* cm~^. The contours correspond to column densities 5, 10.0, 50.0, 100.0 X 10^*cm~^. Panels (c) and (d) show 
the physical density (S = p/pw) and temperature within a central slice through the HVC in panel (b). 



4.1. /. Dark matter dominated clouds 

Figure 2 shows snapshots of the dark matter dominated 
HVC after 3 x 10® years which is approximately 20 crossing 
times {tcr ~ 2R/vw) at a wind velocity Vw = 200 km/s. 
The physical size of the box is 12 kpcx 12 kpcx24 kpc, and 
the number of cells used 128 x 128 x 256, giving a nominal 
numerical resolution of 90 pc. We ran several test simu- 
lations with 256 x 256 x 512 cells and found no difference 
with the lower resolutions runs. The grey scale shows the 
projected surface density of neutral gas and we plot col- 
umn density contours of cold HI. We define a tail as visible 
if it has column density higher than lO^^cm"^. For wind 
densities < 2 x 10~^cm~^ we find no visible tails, although 
very weak features with Nhi ^ 10^® cm~^ are present. 

Panel (b) of Figure 2 shows the same dark matter dom- 
inated HVC moving though a medium an order of magni- 
tude denser. Although this may be unrealistic for a Galac- 
tic halo component, it does give rise to visible tails of gas. 
All of the neutral gas has been pushed away from the dark 
matter potential, as expected from the application of the 
Gunn and Gott stripping criteria. This material will even- 
tually be completely stripped away and evaporate within 
the hot component. 

In all cases that we simulate we find a similar behaviour. 
A weak bow shock is visible in front of the cloud that pen- 
etrates the IGM and increases its local temperature. The 
neutral gas at the front edge of the HVC is compressed 
and heated, whilst a region of low pressure is created im- 
mediately behind the cloud. Panels (c) and (d) of Figure 
2 shows the rich structure and detail that can be cap- 
tured with our code. If the wind is not strong enough to 
cause ram-pressure stripping, viscous/turbulent stripping 
will continue to remove gas, however, the timescale for this 
process is fairly long, hence the tails can be quite weak. 



4.2. //; Pure gas clouds 

Figure 3 shows the evolution of the pure gas clouds. In 
each case, the physical size of the box is 180 pcxl80 pcx 
360 pc the number of grid cells used is 128 x 128 x 256, 
giving a numerical resolution of 3 pc. The evolution and 
features are similar as the dark matter dominated clouds; 
when the wind is dense enough then visible tails of stripped 
HI are created. 

If the initial column density is too low, the pure gas 
clouds are compressed at the leading edges and the veloc- 
ity of the entire cloud decreases as it moves towards its 
terminal velocity. The main difference between the HVCs 
in Figure 3(a) and 3(d) are their initial gas densities, which 
is an order of magnitude higher in Figure 3(d). The IGM 
density and velocities are the same, but the denser HVC 
shows visible tails whereas the tails in Figure 3(a) are be- 
low 10^^ cm^^ and would not be observable. 

The mass ratio of gas in the head and tails is typically 
around 10:1. This result appears very robust regardless of 
the parameters used in the simulation and is lower than 
the 3:1 ratio quoted by Briins et al. (2000). The charac- 
teristic lifetimes of the pure gas clouds are quite short and 
once they enter a diffuse environment that is sufficiently 
dense to create visible tails, they are rapidly destroyed on 
a timescale of w 10'' Myrs. Therefore the high fraction of 
observed clouds with tails suggests that we are observing 
them soon after they are impacting into a region of high 
density gas. The tails trail directly behind the motion of 
the clouds therefore we might expect an alignment effect 
as the clouds fall directly onto the disk. 

5. SUMMARY 

We have performed hydro-dynamical simulations of 
HVCs moving through a diffuse medium in an attempt to 
explain observations of head-tail structures and compres- 



17.0 17.6 18.3 18.9 19.6 17.0 17.6 18.3 18.9 19.6 17.0 17.6 18.2 18.7 19.3 17.0 17.8 18.7 19.5 20.4 

log(Nn) log(Nn) log(Nn) log(Nn) 

Fig. 3. — Column density maps {N[j[) for several of the pure gas HVCs after 3 Myrs of evolution. The temperature and density of the 
external medium are indicated in each panel, and in each case the wind velocity is 200 km/s except panel (b) which is 400 km/s. The initial 
peak column density in HVCs (a)-(c) is 2 X 10^^ cm~^ whereas model (d) was increased to 2 X 10^" cm~^. The contours correspond to column 
densities 5, 10.0, 50.0, 100.0 x lo'**cm-2. 



sion fronts. We find that these features are reproduced 
in either pure gas clouds or gas embedded within a dark 
matter halo, as long as the wind density is higher than 
10~^ cm""^. This is supported by a simple application of 
Gunn & Gott's (f 972) stripping criteria. 

A variety of direct and indirect observations suggest that 
the density of gas in the outer halo is n.uj (50kpc)<5 x 10~^ 
cm"-^. Although observational constraints on the halo gas 
density at 100-400 kpc are highly uncertain, our naive ex- 
pectation is that it should decrease at larger radii. The 
ionised gas close to the Galactic disk can reach much 
higher densities, 10~^ cm~'^, therefore a Galactic foun- 
tain mechanism can reproduce these observations. If the 
halo gas density increases between 50 and 300 kpc by a 
factor of two then the observed head-tail structures could 
also be associated with infalling extra-galactic dark matter 
dominated clouds. 

The lifetime of the neutral tails in our simulations is 
« 10^ years after which time the column density drops 
below the observational threshold. Since the observations 



suggest that approximately 20% of HVCs show evidence 
of tails or compression fronts, a timing problem may ex- 
ist since we must be witnessing the destruction of a large 
fraction of the HVC population at the current epoch. The 
Galactic fountain can circumvent this problem since clouds 
are continuously created and destroyed. The high fraction 
of clouds with signs of interaction is expected since clouds 
are only visible as they condense and fall back onto the 
disk at which point they are susceptible to stripping. 

Our simulations show that all HVCs should show tails 
of stripped gas at lower column densities <10^^ cm~^ 
which may be detectable with higher resolution observa- 
tions. One may also detect Ha emission from the sec- 
ondary shocks that heat the front of the cloud. 
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